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1 Introduction 

^ I Recently, a lot of effort has been paid to the efhcient computation of Kriging predictors 

when observations are assimilated sequentially. In particular, Kriging update formulae 
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enabling significant computational savings were derived in ^|3l|2]. Taking advantage 
of the previous Kriging mean and variance calculations helps avoiding a costly (n + 1) x 
(n + 1) matrix inversion when adding one observation to the n already available ones. In 
addition to traditional update formulae taking into account a single new observation, |2] 
also proposed formulae for the batch-sequential case, i.e. when r > 1 new observations 
are simultaneously assimilated. However, the Kriging variance and covariance formulae 
given without proof in [2] for the batch-sequential case are not correct. In this paper 
^ i we fix this issue and establish corrected expressions for updated Kriging variances and 

Cn ' covariances when assimilating several observations in parallel. 
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2 Kriging update formulae for the parallel case 
2.1 notations and formulae 



In this section we give a counter-example for the Kriging variance update formula pro- 
posed in [2]. Let us first introduce the notations of this paper and recall these formulae. 
^ , When n observations, at locations xi, . . . ,x„ € D, of a real valued random field Z in- 

dexed hy D C W^ are available, [2] denotes respectively by Z{-x.)n data and o"^(x) the 
Kriging mean and variance at x S D based on these n observations. The corresponding 
Kriging covariance is denoted by (T„ : (x, y) S D^ -^ crn{x, y). If /c > 1 new observations, 
at points x„+i, . . . , x„+fc, are available, the parallel Kriging update formulae given in [2] 
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follow: 

k 
Z{'X.)n+k data = Z{-X.)n data + ^^^^ A„_|_j|„+fc(x) f Z(x„+j) — Z{^n+i)n data) , (1) 

k 

(rl+ki^) =o-^(x) - ^A^+i|„+^(x)cj^(x„+i) , (2) 

i=l 

o-n+fc(x,y) =cr„(x,y) - ^ A„+j|„+fc(x)A„+i|„+;t(y)o-.^(x„+j) (3) 

where A„+j|„_|_fc(x) denotes the Kriging weight of Z(x„-|_j) when predicting Z{x.) relying 
on n + A: observations. In [2], these formulae are proven only for k = 1. For k > 1, 
they are quickly justified with an argument based on Pythagoras' theorem. In fact, the 
assumptions for using this theorem are not satisfied, as we will detail here. We now give 
a counter-example for Equation ([2]) and then propose fixed expressions. 

2.2 A counterexample for Equation ([2]) 

Let us consider for the sake of simplicity a case where Z is a one-dimensional centered 
random process indexed by -D = [0,1], with covariance kernel C{x,y) := cro(x,y) = 
inm{x,y) (Wiener process, or Brownian Motion). Note that even though this process 
is non-stationary, Simple Kriging equations are applicable to it (Written in the general 
case of a non-stationary covariance kernel, as allowed in [2]'s results and notations). 
Assuming n = initial observations and k = 2 new observations at the points xi = ^, 
X2 = 1, the Kriging weights for a prediction at x = | write: 



4(x) = C(x,x) - (A,p(i),A2p(rc)) (';.t]'?} I^i!!'!'!) (A,p(i), A2|2(i)f = ^ 



leading to a Kriging mean Z(x)„+fc data = \Z{xi) + ^Z{x2), and to a Kriging variance: 

Ha 

\C{x2,Xi) C{X2,X2)^ 

Now, using Eq. ([2]) would lead to 

o"i(a:) = o-g(x) - A^|2(x)o-g(xi) - Xl^2ix)'^oi^2) 
_ 3 11 1 _ 3 
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and we see that the two expressions do not match. More precisely, the conditional 
variance obtained by using Eq. ([2]) is too large, as if the contribution of xi and X2 to 
decreasing the variance would have been underestimated. As we will see now, this is 
indeed the case, and the missing part is related to conditional covariances. 



2.3 Corrected Kriging variance and covariance update formulae 

In this section, we establish corrected update formulae that may be used in the batch- 
sequential case instead of the formulae ([2]) and ([3]). To improve the readability of the 
properties and their proofs, we adopt the following simplified notations: 

• Xoid := {xi, . . . ,x„}, and Xnew := {xn+i, • • • ,x„+fc}, 

• Zoid := {Z{xi), ..., Z(x„)), and Znew := (^(x„+i), . . . , Z(x„+fc)), 

• '^ncw.oldlx) := (Ai|„+fc(x), . . . , A„|„_|_jt(x)) , 

• -^ncw.ncwix) := (A„^_i|„4.fc(x), . . . , A„_|_fc|„_|_fc(x)) , 

• (T^jj(x) := (T^(x), cr^j;„(x) := (T^^^(x), and similarly for the conditional covariances. 

• The conditional covariance matrix of Znew knowing Zqm is denoted by Snew; 

For conciseness and coherence, Z{x.)n data and Z(x)„_|_fc data are also denoted by Z(x)oid 
and Z(x)new, respectively. The corrected update formulae are given below: 

Proposition 1. (Corrected Kriging update equations for the parallel case) 

■^yXjnetu — ^x^jold ~r '^new,new\^j y^new ^\-^new)old) \^) 

'^netuyX, yj — (Jold\^Ty) ^new,new\^) ^new^new,new\y ) (Oj 

The proofs above use a Gaussian assumption on the field Z, enabling a convenient 
interpretation of the Kriging mean and variance in terms of conditional expectation and 
variance. A key property is that even though the conditional distribution interpretation 
doesn't hold in non-Gaussian cases, the formulae for the Kriging weights remain valid 
whatever the assumed square-integrable distribution for the field; we are just using here 
that best linear prediction and conditioning coincide in the Gaussian case. Furthermore, 
the results hold as well for the cases of Ordinary and Universal Kriging (written in 
terms of covariances, for a square-integrable field, not with variograms), relying on their 
well-known Bayesian construction with an improper prior on the trend coefficient (s). 

Proof. The first equation follows from an application of the law of total expectation: 



Z(x)oid = E[Z(x)|Zoid] 

= E[E[Z(x)|Zold,Znew]lZold] 



E[ Anew, old (x) ZqM + Anew, new (x) ZnewJ|Zo 



'^new,oldiXj ZjQjd + Anew, new (Xj -^v-^ncwj 
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Anew,oldvXj AqM i Anew, new (Xj ^new Anew,newvXj Zjnew i Anew, new (,Xj ^ (.-'^newjold 
^ (Xjnew Anew,new(,Xj (^Zjnew ^ (.-'^newjoldj 



Similarly, using the law of total variance delivers: 

C72id(x) = Far [Z(x) I Zoid] 

= E[yar[Z(x)|Zoid,Znew]lZoid] + V^ar[E[Z(x)|Zoid,Znew]|Zoid] 

= Var[Z {x)\Zold, Znew] + ^a?'[Ancw,old(x)^Zold + Anew,ncw(x)^Znew|Zold] 



'-''nowV-Xj ~r '^ncw,ncwvXj ^ 



now'*ncw,ncw i, 



which proves the second equation. The third equation comes with the same method, 
using the law of total covariance. D 

Proposition 2. (Kriging weights expressed in terms of conditional covariance) 
Corollary 1. (Corrected Kriging update equations in terms of conditional covariance) 

Z(pi.)new = ^(x)oM + 0'oM(X„ew5 x) S„g^(Z„ew ~ Z[X.new)old) (8) 

^lewi^) = (^lldi^) - (^oldO^new, :>i.)'^T.~}^aold{'^new, x) (9) 

0"nej«(x,y) = 0-oid(x, y) -0-oM(X„eM, x) S~g„(ToM(X„e«„ Y) (10) 

Proof. Using the orthogonal projection interpretation of the conditional expectation, 



Z(x) = E(Z(x)|Zold, Znew) + ^(x) - E(Z(x)|Zold, Z^cw) 
— '^new,old (,Xj ^old ~r Anew, new (Xj ^new ~r £ , 

with £ centered, and independent of Zqm and Znew- Let us now calculate the conditional 
covariance between Z{x.) and Znew knowing the observations Zqm: 

0-old(Xnew,x) := COW ( Znew, -^(x) |Zold) 

= + GOV ( Znew, Anew,new(x) Znew | ZqM j + COt;( Znew, e|Zold) 
— ^new Anew, new (Xj ~r COV[^£, Zinewl^oldj 

Noting that cou ( Znew, e I Zqm) = 0, the latter equation proves Proposition [2j Eqs. [HI El 
and M of the corollary directly follow by plugging in Eq. [2] into Eqs. [H El O D 

2.4 Interpretation on the counter-example 

The main difference between Equations [2] and [5] is that Eq. [5] takes into account condi- 
tional Kriging covariances. Coming back to the counter example, Eq. [5] delivers: 

(^Ux) = CToix) - Ai|2(2:)cro(xi) - A2|2(a;)cro(2:2)-2Ai|2(x)A2|2(x)cro(xi,X2) 

_ 3 11 1 11 _ 1 

~4~42~4~ 42 ~8 



which is the correct variance, as obtained earher with regular Simple Kriging equations. 
Note that Eq. [2] would be correct if the matrix Snew were diagonal, which has no par- 
ticular reason to happen in practical applications. The Pytagorean theorem invoked in 
[2], enabling to decompose a variance into a sum of (weighted) variances, cannot hence 
be applied in the general case because of non- independence between Kriging residuals. 

3 Conclusion 

In this paper we derived corrected Kriging update formulae, fixing the incorrect formu- 
lae [2] and [3] given in [2] , where a part of the update involving conditional covariances 
had been neglected. A simple interpretation of the new update formulae is that they 
correspond to Simple Kriging equations, where the underling covariance kernel would 
be the conditional covariance kernel before update, i.e (Toid(-, •)• These formulae enable 
important computational savings (avoiding a large matrix inversion) and may perfectly 
be adapted to further frameworks such as co-Kriging. 
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